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Abstract 

In an extensive numerical investigation of nonintegrable translational motion of discrete breathers 
in nonlinear Schrodinger lattices, we have used a regularized Newton algorithm to continue these 
solutions from the limit of the integrable Ablowitz-Ladik lattice. These solutions are shown to be 
a superposition of a localized moving core and an excited extended state (background) to which 
the localized moving pulse is spatially asymptotic. The background is a linear combination of 
small amplitude nonlinear resonant plane waves and it plays an essential role in the energy balance 
governing the translational motion of the localized core. Perturbative collective variable theory 
predictions are critically analyzed in the light of the numerical results. 
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Discrete breathers are spatially localized, time periodic solutions of homoge- 
neous nonlinear lattices, which have been recently observed in experiments on a 
variety of physical systems (magnetic solids, arrays of Josephson junctions, cou- 
pled optical waveguides and photonic crystals). Though many of the properties 
of discrete breathers are today well characterized, the question on their mobility 
remains under controversy, due to the radiative losses unavoidably associated to 
the traslational motion of the localized pulse in generic (nonintegrable) systems. 
We address here this problem in an important class of nonlinear lattices: the 
discretizations of the Nonlinear Schrddinger Equation. Our results show that 
exact mobile breather solutions ride over an extended excited state of the lat- 
tice, which we fully characterize. Moreover this background p lays an essential 
role in the energy balance required for exact nonintegrable mobility. 



I. INTRODUCTION 



Nonlinear lattices have become the subject of a considerable multidisciplinary interest, 
with applications in physics subdisciplines as diverse as biophysics (myelinated nerve fibers 
Q, DNA 

biopolymer chains |3[), nonlinear optical devices (photonic crystals jj| and 



waveguides 



Bose-Einstein 



and Josephson effect Q| (superconducting devices 0,0], 
condensates |lflL IllL Il2j|). among others. From a theoretical perspective they have been 
progressively recognized not as mere discretizations (unavoidable for numerical computa- 
tions) of nonlinear continuum field equations, but as a target of interest in their own right, 
due to the distinctive features associated with discreteness, whose relevance to experimental 
features have been largely established. 

More specifically, among the variety of observable nonlinear behaviors of lattice dynamics, 
the phenomenon of (discrete breathers) nonlinear localization in lattices 3| has received 
attention in both experimental and theoretical research during the last several years. (Non- 
topological) discrete breathers are exact spatially localized, time-periodic solutions. Due 
to discreteness the plane wave spectra are bounded, thus making posible the absence of 
multi-harmonic resonances of the exact discrete breather solution with extended modes. 
The combination of nonlinearity and discreteness is sufficient for the physical existence of 
discrete breathers, DB's for short, resulting in its generality and broad interest. The reader 
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may find in Q| a recent multidisciplinar survey of current research on the subject. 

Our primary concern here is the issue of DB's mobility. The translational motion of 
discrete breathers introduces a new time scale (the inverse velocity), so generically a moving 
breather excites resonances with plane wave band spectra. This fact poses no problem to the 
persistence of localization when the lattice dynamics is governed by power balance (forced 
and damped lattices the emitted power is exactly compensated by the input from the 

homogeneous external force field, during stationary breather motion. However, for generic 
(nonintegrable) hamiltonian lattices, the radiative losses would tend to delocalize energy, 
and some energy compensating mechanism is needed in order to sustain exact stationary 
states of breather translational motion. From the (particle) perspective of collective variables 
theory, the localized breather experiences a periodic Peierls-Nabarro potential function of 
its position, so that the motion of the localized field oscillation over this landscape should be 
expected to induce the emission of radiation at the expense of translational (and/or internal) 
breather kinetic energy, which thus would unavoidably decay on time. 

To address the problem, a reasonable strategy is to use precise numerical techniques on 
adequately general models, with the hope that they may pave the way to further physical 
(and mathematical) insights. Our chosen model belongs to an important class of nonlinear 
lattice models: the discretizations of the continuum nonlinear Schrodinger equation, here 
referred to as Nonlinear Schrodinger (NLS) lattices. First of all, they are ubiquitous in 
models for polaronic effects in condensed matter, nonlinear optical technologies, and the 
physics of Bose condensate lattices and superconducting devices, where nonlinear localiza- 
tion is currently studied. Second (a very convenient technical advantage), this class contains 
an integrable limit (the Ablowitz-Ladik lattice, see below) having exact moving discrete 
breathers, wherefrom perturbation (collective variable) theories have been developed in sup- 
port of exact (or very approximate, in the least) nonintegrable mobile DB's. This allows a 
detailed theoretical analysis of the numerical results, as well as eventual feedbacks for useful 
(and currently used) theoretical concepts and perspectives. In the subsection II Al of this 
introduction we present the (three-parameter) NLS lattice that we have studied, originally 
introduced by Salerno 16^ . 

The numerical techniques employed are summarily introduced in section II. We stress here 
the unbiased character of this numerical procedure which, unlike other techniques, is not 
based on ansatze on the expected functional form of the exact solution sought. In essence, 
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the procedure uses a regularized Newton continuation method for operator fixed points, and 
it only requires a good starting set of parameter values where the exact solution is known. 
In our case this is provided by the integrable Ablowitz-Ladik limit of the NLS lattice, from 
which adiabatic continuation of the two-parameter (core frequency a>& and velocity v 5) family 
of moving Schrodinger breathers is performed. 

The main numerical facts are shown in section III. The numerical solutions are found to be 
(up to numerical precision) the superposition of a traveling exponentially localized oscillation 
(the core), and an extended background, which is a linear superposition of finite amplitude 
nonlinear plane waves Aexp[i(kn — cut)] (see Fig. [TJ. The amplitudes of these resonant 
nonlinear plane waves are observed to differ typically by orders of magnitude, so that only 
a small number of them are relevant for most practical purposes. They fit well simple 
theoretical (thermodynamic limit) predictions based on discrete symmetries requirements. 
Contrary to the exact immobile breather solution (space-homoclinic and time-periodic orbit), 
which asymptotically connects the rest state (vacuum or ground state) of the lattice with 
itself, each exact mobile localized solution is instead homoclinic to a specific lattice state of 
extended radiation. In other words, exact stationary mobility of discrete breathers requires 
an extended excited state of the lattice. A preliminary account of some of the numerical 
results of this section were reported in |l7| . 




FIG. 1: Time evolution of \& n \ 2 profile of a mobile discrete Schrodinger breather. The frequency 
of the solution is uif, = 5.050 and the velocity is vj, = 0.804. Note that the background is composed 
by a single plane wave with amplitude A. The nonintegrable parameter of eq. ((7)) is v = 0.2. 

In section IV we analyze the numerical results in the light of collective variable the- 
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ories, correlating them with the main theoretical predictions of this successful (however 
incomplete) physical perspective. In particular, the existence of Peierls-Nabarro barriers to 
translational core motion is confirmed, and its subtle relation to the background amplitude 
is discussed. We present also numerical confirmation of the existence of exact oscillating 
anchored breathers, whose extended background is much smaller than those of traveling dis- 
crete breathers of the same internal frequency Ub- Along with the discussion in this section, 
a physical interpretation of the role of the interaction background-core in the energy balance 
emerges, paving the way to a satisfactory integration of the results into a collective variable 
theory. 

Finally, in section V, after summarizing the main conclusions drawn on discrete 
Schrodinger breather mobility, we briefly trace some interesting open questions for fur- 
ther research, notably the approach to irrational breather time scales ratios, the study of 
multibreather solutions (two-breather collision processes, trains of moving breathers), and 
the coupling to both thermal and nonthermal (e.g. elastic) degrees of freedom, where the 
numerical tools and results presented here can find further applications. 



A. NLS lattices 

The standard discrete nonlinear Schrodinger (DNLS) equation is the simplest dis- 

cretization of the one-dimensional continuous Schrodinger equation with cubic nonlinearity 
in the interaction term, i.e., 

i$ n = -C($ n+1 +$„_l) -7|$n| 2 $n • (1) 

In this expression $> n (t) is a complex probability amplitude, the parameter C amounts 
the nearest neighbor coupling, and 7 is the strength of the nonlinearity. The self-focussing 
effect of local nonlinearity balanced by the opposite effect of the dispersive coupling makes 
possible the existence of localized periodic solutions (breathers) of the discrete field, where 
the profile of |$ n | decays exponentially away from the localization center: 

= |$J esxp[kj b {t))] . (2) 

In the uncoupled limit C — > of the DNLS equation, also known as the anti-integrable or 
anti-continuous limit, discrete breathers can be easily constructed by selecting a periodic 

5 



oscillation § no (t) of frequency u>b at site no and $ n = for n ^ n . These solutions can be 
uniquely continued to nonzero values of the coupling C, and constitute the one-parameter 
family of immobile on-site breathers of the DNLS equation. 

Unfortunately the continuation from the uncoupled limit does not provide solutions where 
the localization center moves along the lattice with velocity v^, i.e, mobile discrete breathers. 
On the other hand, there is an integrable lattice as a limit of the nonlinear Schrodinger class 
that posseses this type of mobile solution. That is the one discovered by Ablowitz and Ladik 
in [ljj: 

i$ n = -C($ n+1 + $ n _0 [l + ||$ n | 2 ] , (3) 

where, again, C and 7 account for the strength of the coupling and the nonlinearity, respec- 
tively. The integrable Ablowitz-Ladik equation, A-L for short, possesses a two-parameter 
family of exact moving breather solutions of the form 

$ n (t) = y^-sinh/3 sech[i3(n - x Q (t))] x 

exp [i(a(n - x (t)) + Q(t))] . (4) 

The two parameters of this breather family can be chosen to be the breather frequency u>b 
and velocity t>&, 

2 sinh (3 sin a , . 

v b = x = (5) 

u>b = Ct = 2 cosh (5 cos a + avb , (6) 

where — n < a < tt and < (3 < 00. The A-L moving breather (instantaneous) profile 
interpolates between the rest state $ n = of the lattice (at n — > ±00) in an exponentially 
localized region around Xo(t), while traveling with velocity Vb- 

The connection between the integrable (though physically limited) A-L equation and 
the physically relevant (though nonintegrable) DNLS equation is provided by the model 
originally introduced by Salerno in 

i$ n = -($ n+1 + $„_!) [C + /i|$ n | 2 ] - 2z/$ n |$ n | 2 . (7) 

This lattice provides a Hamiltonian interpolation between the standard DNLS equation (jlj). 
for u = and v = 7/2, and the integrable A-L lattice when « = 7/2 and v = 0. The 
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Hamiltonian of the Salerno equation is given by 



n 



- C^($ n $„ +1 + $ n $ n+1 )-2-^|$, 

n n 

+ 24^1n(l + /i|$„| 2 ), 



(8) 



which contains the A-L and DNLS Hamiltonian for the above limits. In addition to the 
Hamiltonian, this equation possesses, for any value of the parameters, the following con- 
served norm: 



In the following we will set the value of 7 = 2 (as usual) and consider the coupling strenght 
C = 1 in eq. (JZJ). 

The continuation of the family of mobile discrete breathers from the A-L integrable limit 
allows numerical observations of the interplay between the integrable term, weighted by the 
parameter /1, and the nonintegrability, weigthed by v. 

II. DISCRETE BREATHERS NUMERICS 

We introduce here the numerical techniques that we have used. As a whole, one could 
refer to them as the (SVD-) regularized Newton method. They do not constitute a novel 
method in "discrete Breather numerics", as they have been already used, e.g. in j^] to 
refine movin g b reathers of Klein-Gordon lattices obtained by other numerical means (see, 



by contrast, [21]). From the methodological side, what is novel here is the systematic use of 



them in the investigation of the family of moving Schrodinger breathers reported below in 



To some extent, the presentation here is self-contained. First, in III Al we introduce the 
notion of (p, q) resonant solution, providing some illustrative examples. The (SVD) regular- 
ized Newton algorithm is presented in III Bl and finally in III CI we briefly explain the basics 
of Floquet stability analysis. 

A. Discrete space-time symmetries: (p, q) resonant states 

If a frequency u b = 2n/T b is given, we will say that a solution $ = {$„(£)} is (p, q) 
resonant with respect to the reference frequency u b , if the following condition holds, for all 




(9) 




EH 
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n and t: 

$ n (t) = $ n+p (t + gT b ) . (10) 

After q TVperiods, these solutions repeat the same profile but displaced by p lattice sites. 
In more technical terms, these (p, q) resonant solutions are fixed points $ of the operator 

C p T q = M (11) 
(M — J) $ = , (12) 

where £ and T are, respectively, the lattice translation and the T^-time evolution operator 

£{$ n (t)} = {$„ +1 (t)} (13) 
T{$„(t)} = {<£> n (t + T b )} . (14) 

We now consider some examples of (p, q) resonant solutions with respect to the frequency 
cj&; the first example is simply provided by the family of plane wave solutions of eq. (J7J): 

$ n (t) = Aexp[i(kn - cot)} . (15) 

It is easily seen, by inserting (|T5|) in eq. ( |7j) . that the values of a;, and |A| define a surface 
in the three-dimensional space, the nonlinear dispersion relation surface u(k,A) (see Fig. 
0): 

u = -2[1 + fi\A\ 2 } cos k -2v\A\ 2 . (16) 

Note that due to the nonlinear character of the eq. (J7|), the frequency uj depends on both 
wave number k and amplitude \A\ of the plane wave. 

One can easily determine those plane waves that are (p, q) resonant with respect to u^. 
the eq. (jl()j) imposes the following condition on uj and k 

to 1 / p 



, -k-m , (17) 

where m is any arbitrary integer. These planes in the 3-d space (uj, \A\, k) intersect the 
dispersion relation surface at (in general) several one-parameter families (branches) 
in the first Brillouin zone (—it < k < ir). 

If we are not interested in unreasonably large (and not interesting) amplitude values \A\ 
of the plane waves, the number of branches is finite: one can see that for fixed values of all 
the parameters (p, q, u>b, u, fi), there is a finite number of branches in the limit \A\ — ► 0; 
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FIG. 2: Plot of the nonlinear dispersion relation surface of nonlinear plane waves, eq. (|16j) . as a 
function of the amplitude A and the wave number k of the plane wave. The values of \x and v are 
fixed to 0.5. 

there is also a well defined (parameter dependent) treshold value of the amplitude at which a 
pair of new branches (tangent bifurcation) appear (i.e. these plane waves can only resonate 
with u>b for amplitudes above some threshold value). 

Thus, by a suitable bounding of the amplitude, for each couple (p, q) one finds a finite 
number, s, of branches of (p, q) resonant plane waves. (Note also that this number diverges 
when p/q tends to an irrational). 

A different, and highly nontrivial, example of (p, q) resonant solutions is provided by the 
solitary waves (jH) of the A-L lattice. From eq. © it is clear that the choice litv^fiQi = p/q 
selects a (p, q) resonant solitary wave with respect to the frequency ujb- The set of velocity 
values of resonant A-L breathers is dense and any A-L moving breather is a limit of some 
sequence of resonant ones. Note also that immobile breathers are (0, 1) resonant with respect 
to the frequency u b . 

In the integrable limit, the plane waves and the A-L breathers are both exact independent 
solutions. Integrability makes possible that the initial localization of energy is maintained 
with time evolution, without decaying away by exciting radiation. It is a well established 
result that (even far away from this integrable limit) immobile discrete breathers remain 
exact solutions of the lattice dynamics. Our concern in the next sections is the question of 
moving discrete breathers away from integrability in eq. (J7J). In order to study them, we 
will focus on (p, q) resonant solutions. The motivating of this restriction comes from its 
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accessibility to numerics. First we will motivate the numerical (Newton) method that allow 
us to study these solutions with an adequately high precision. 



B. Newton continuation 



A well-known numerical procedure to o 
is the Newton continuation 

0, Q, Q 



)btain exact periodic solutions of nonlinear lattices 
241 ] . The different practical implementations of this 
procedure work very successfully when, for example, one obtains numerically exact immobile 
discrete breathers of eq. (JJJ), from the uncoupled limit fi = and C — 0, where exact periodic 
discrete breathers are trivially constructed. 

The iteration of the Newton operator T converges rapidly to its fixed point (i.e. the 
solution to be computed) provided the starting point, $°, is close enough, and the solution 
of the following system of linear equations is a well-posed problem: 

(DT — l)($ n — $ n+1 ) = [T — J]$ n , (18) 

where DT is the jacobian matrix of the Newton operator, and $ n (the ra-th iteration solution 
of ()18|)) converges quadratically to the fixed point solution. By adiabatic change of a model 
parameter, one constructs a uniquely continued exact fixed point solution for each parameter 
value, using each time, as starting point of the Newton iteration, the solution previously 
computed. 

The matrix (DT —1) must be invertible, in order to uniquely compute $ n+1 . Degeneracies 
associated with the +1 eigenvalues of DT, if any, have to be removed in order to obtain a 
unique fixed point solution. When continuing immobile (periodic) discrete breathers of eq. 
(J7J), a convenient prescription is commonly used, namely to restrict the operator action to 
the subspace of time-reversible solutions 0,0]. This provides a practical way of removing 
degeneracies, allowing unique continuation of immobile discrete breathers. 

For the continuation of (p, q) resonant solutions (of which periodic solutions are only the 
particular case p = and q — 1), one has to use M. = C p T q as the Newton operator. One 
has also to deal with the degeneracies of A4, and imposing time-reversibility could, in this 
case, be too restrictive. 

A well-known solution to the problem of removing degeneracies when no clear restrictions 



are available, is provided by the so-called singular value decomposition (SVD) 



20 



22 
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of the matrix (DC p T q - 1) 



(DC p T q - 1) = J = PVQ 



(19) 



where P, V and Q are 2N x 2N square matrices. P and Q are orthogonal matrices and V 
is diagonal (vjSij) with possibly null (zero) elements, called singular values, associated with 
the null space of J (the subspace that is mapped to zero Jx — 0). The columns of P whose 
same-numbered elements Vj are nonzero are an orthonormal set of basis vectors that span 
the range of J (the subspace reached by this matrix). The rows of Q whose same-numbered 
elements Vj are zero are an orthonormal basis for the null space of J. One can numerically 
use this SVD decomposition, checking the (numerical) vectors spanning the null space to 
identify degeneracies, and using at iteration steps the pseudoinverse matrix 



where V~ l is diagonal with elements 1/vj for Vj and for Vj = 0. 

As a judicious test of our numerical codes, we have used both procedures (reduction to 
time-reversible subspace and SVD decomposition) to obtain immobile discrete breathers of 
the Salerno model. Both agree, up to the highest possible accuracy, from the uncoupled 
limit (one-site and two-site breathers) up to the A-L limit (and viceversa). 

This test serves also to provide further confirmation of an important and well-known 
theoretical result. At the integrable A-L lattice, one-site and two-site immobile breathers 
are but two particular choices of the continuous one-parameter (xq, the localization center) 
family of immobile solitary waves, i.e. constant Xo(t) = n or n + 1/2 respectively, in eq. (jljl. 
The well-known result, confirmed by our numerics, is that away from the A-L limit only these 
(one-site and two-site) immobile discrete breathers persist under adiabatic continuation. No 
immobile breather centered in between exists. For positive values of the parameter u, the 
one-site immobile one has a lower value of energy 7i, and it is a linearly stable solution, while 
the energy of the two-site breather is higher and it is linearly unstable. The relative situation 
is reversed for negative values of v. This result can be interpreted as the emergence of a 
(Peierls-Nabarro) potential function of the breather center Xq, which destroys the continuous 
degeneracy of immobile breathers, leaving only two of them per lattice unit, namely those 
localized at maxima and minima of the Peierls potential. This interpretation, which is 
captured in the theoretical framework of collective variable approaches, turns out to play 



Q*V- l P* , 



(20) 
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a central role in building up the physical interpretation of the numerical results on mobile 
discrete Schrodinger breathers, below in section ITTT1 

The numerical integration of the equations was performed using a 4 th order Runge-Kutta 
scheme with time step T^/500. The convergence criterion for the fixed point solution is 
that the value of ~~ ^] < ^ >n+1 )jl i s l ess than N ■ 10~ 16 (where N is the size of the 

lattice). The typical size of the lattices was taken between 100 and 200 sites depending on 
the characteristics of the solution considered, as we will explain in section IIHI 



C. Floquet stability analysis 

A very useful outcome of the numerical Newton method of computing solutions of eq. (J7J) 
is the jacobian matrix of the Newton operator, usually called the Floquet matrix F. This 



matrix is the linear operator associated with the linear stability problem 27[ of the fixed 
point solution. 

Indeed, the jacobian F of the Newton operator M. 

F = DM (21) 

maps vectors in the tangent space of the solution (small initial perturbations e*(0) of the 
fixed point solution) into their T^-evolved vectors, i.e. 6(Tm), after a period of M.. That 
is: 

e(T M ) = jFe(0) , (22) 

The Floquet matrix of a Hamiltonian system is real and symplectic, so the Floquet 
eigenvalues A come in quadruplets, A, 1/A, A, 1/A. The necessary condition for the stability 
of the solution is that all the eigenvalues lie on the unit circle of the complex plane, |A| = 1. 

To illustrate the Floquet analysis of (p, q) resonant solutions of the NLS lattice (J7|), we 
now obtain the Floquet spectrum of modulational instabilities of a (p, q) resonant plane 
wave, 

$„(t) = Aexpi(kn - cut) . (23) 
(The modulational instabilities of plane wave solutions of nonlinear lattices have been ana- 



lyzed in |2i 



2*. 
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One has to investigate the evolution of small perturbations, in both amplitude and phase, 
of the plane wave 

$ n (t) = (A + I n ) exp i(kn -ut + tp n ) , (24) 

where we assume that the perturbation parameters are small compared with those of the 
plane wave solution. Introducing expression (J24j) in (JJJ) and considering the following form 
for the perturbations {I n ,<p n }: 

I n (t) = I exp i(Qn — Qt) 

ip n (t) = ip exp i(Qn — Vtt) (25) 

we obtain the dispersion relation for the perturbation parameter Q: 

[Q - 2(1 + fiA 2 ) sink sin Q} 2 = 16(1 + fiA 2 ) x 

sin 2 Q/2cosA;[(l + /j,A 2 ) sin 2 Q/2 cos k 

- fiA 2 cosk - v A 2 }, (26) 

From the above expression one derives the values of Q(A, Q, k; u, n) for the modulational 
perturbations. When the parameter Q has a nonzero imaginary part, i.e. the right-hand 
side of (J26|) is negative, the plane wave (A, k) becomes unstable under the corresponding 
modulational (Q) perturbation, whose amplitude will grow exponentially fast in the linear 
regime (tangent space). 

Modulational perturbations (|23|) correspond to eigenvectors {I n , <p n } of the Floquet ma- 
trix: 

I n (t + T M ) = exp(-iQT M )I n (t) (27) 
ip n (t + T M ) = exp(-iflT M )ip n (t) (28) 

with associated Floquet eigenvalues exp(— if2T^). The real part of Q gives the angle in the 
complex plane, 

Ofi*, = -K(tyT M , (29) 
while the imaginary part gives the modulus of the Floquet eigenvalue, 

|A| = exp(3f(fi)T JW ) , (30) 
thus providing the information about the linear stability of the solution. 
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The distribution of angles and moduli in the Floquet spectrum of the modulational in- 
stability can be obtained from eq. (|2T)|) by taking the real and imaginary parts of Q: 

&(ft) = 2(1 + fiA 2 ) sin k sin Q (31) 

3(fi) 2 = -16(1 + iiA 2 ) sin 2 Q/2 cos k x 

x [(1 + /jA 2 ) sin 2 Q/2cosA; 

- /iA 2 cosA; - iM 2 ] . (32) 



A*=0.01 
A,=0.02 
A,=0.03 
A,=0.04 
A =0.05 

-0.5 0.5 1 1.5 

e Floq I T M 

FIG. 3: Plot of the modulus of the unstable Floquet eigenvalues |A| (corresponding to the positive 
values of G(O) in eqs. l|3*Tj) and (|32|lh versus the Floquet angle, #FZog- Both quantities are conve- 
niently normalized to the period of the map Tm- The amplitude of the excursion of |A| and the 
range of values of 9 Floq f° r which |A| > 1 grow as the amplitude A of the plane wave is increased. 
The parameters in eq. (J7J) are [i = v = 0.5 and the wave number of the plane wave is k = 0.5. 

In Fig. 01 we represent the modulus of the unstable eigenvalues as a function of the 
Floquet angle for the spectrum of a (p, q) resonant plane wave, taken as an example to 
visualize the non-point-like character of the instability in the Floquet spectrum in the ther- 
modynamic limit. Note that there is no plane wave harmonic instability {9pi oq = 0) due to 
this mechanism of modulational instabilities. 

A numerical computation of the Floquet spectrum of a plane wave (with arbitrary wave 
number) of a lattice of N = 400 sites, with periodic boundary conditions is shown in the 
complex plane representation of Fig. 0J The instability globes, at angles symmetrically 
placed around zero in this figure, nicely fit the theoretical (thermodynamic limit) values 
obtained from eqs. (j3~Tj) and (|32~j) . 
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FIG. 4: Plot of the Floquet spectra of a plane wave with modulational instability (circles) and the 
theoretical prediction (lines) for the distribution of the Floquet eigenvalues in the complex plane 
given by eqs. (|31|) and (|32|) . The amplitude and wave number of the plane wave are A = 0.1 and 
k = 0.1 • 2ir; the nonintegrable parameter value is v = 0.1 and the lattice size is of 400 sites. 

III. MOBILE DISCRETE SCHRODINGER BREATHERS 

In this section, we show the numerical results on mobile discrete Schrodinger breathers in 
the NLS lattice (|7j). These numerics are computed using the tools explained in the previous 
section. The Newton fixed point continuation requires a good initial guess (meaning that 
the starting initial conditions have to be in a small neighborhood of the fixed point). Very 
close to v = 0, the A-L solitary traveling waves (exact solutions at v — 0) provide good 
starting points. After convergence to the fixed point, we increase adiabatically the value of 
the parameter (Au = 10~ 3 ), and start iteration from the previous fixed point. 

An important step in the numerical method used here, is obtaining a basis for the sub- 
space of (tangent space) vectors with Floquet eigenvalue +1. These are associated to those 
degeneracies (symmetries) that one has to eliminate in order to regularize the linear system 
at each (Newton) iteration step when numerically converging to the fixed point solution. 

Away from the A-L limit, it is known (as reported e.g. in j^J) that only two conserved 
quantities remain generically as dynamical invariants, the Hamiltonian (jHJ) and the norm Q . 
They are respectively associated to the continuous time translation and gauge (global phase 
rotation) invariance. Using the notation Ui = 9ft($j) and v; t = S($j), one easily obtains 
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that (Sui(t) = iii(t), Svi(t) = Vi(t)) is the perturbation associated with time translational 
invariance, while (Sui(t) = Vi(t), Svi(t) = —Ui(t)) is the one with gauge invariance. These 
are, consequently, Floquet eigenvectors with associated eigenvalue +1, and we can easily 
check that they coincide with the (two) basis vectors provided generically (i.e. except at 
special bifurcation values of the parameter, see below in 1111 Cj) by the numerical Singular 
Value Decomposition (J2U)l explained in the previous section. 

In subsection IIII Al we summarize our findings on the generic structure of mobile 
Schrodinger discrete breathers. For this, as explained earlier, we have explored particu- 
lar values for the integers (p, q) and performed continuation of (p, q) resonant A-L traveling 
waves. The variation of the main structural characteristics of the fixed points along the 
continuation parameter v is examined in detail in 1111 Bl for both signs of this parameter. 
Of particular interest are the observed drastic changes in the structure for v ~ —0.3 and 
v ~ —0.39. Then, in IIII CI we show the main conclusions on the stability analysis of the 
mobile Schrodinger discrete breathers, in a sector of the breather parameter space. 

A. The structure of the solution 

In Fig. El we plot the spatial profile of a (1, 1) mobile Schrodinger discrete breather for 
nonintegrability parameter value v = 1.0, and = 2.678. 

A quick inspection of this figure provides a first glance of the general structure of the com- 
puted (p, q) resonant solutions: The fixed point $ is the superposition of an (exponentially) 
localized oscillation (the core) moving on top of an extended background. 

$ = $ core + $b a ckg ■ ( 33 ) 

In other terms, far away from the core localization site n , the solution does not tend to the 
rest state = 0, but to an extended excited state of the lattice , i.e. for \n — n \ 3> 1 

*n(t) = (% ackg )n(t) ± . (34) 

One easily realizes (for example, consider a site very far from no) that the background has 
to be itself (p, q) resonant. This can be quickly checked in our numerics: Indeed, the power 
spectrum S(u) = \ 9ft[<I> n (t)] exp(iu;t)<it| 2 at a site n far from n reveals a finite number 
of s peaks Uj, j = 0, ...,s — 1; one can check that each ujj numerically fits to a branch of 
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FIG. 5: Instantaneous profile of a (1, 1) resonant breather with ajf, = 2.678 and Vj, = 0.426; the 
nonintegrable parameter is v = 1.0 (standard DNLS equation), (a) Real part, (b) imaginary 
part, (c) modulus and (d) phase. The resonant condition for the harmonic composition of the 
background gives the contribution of three plane waves. The existence of these plane waves is 
revealed by the modulation of the extended tail in the modulus profile (c). 

(p, q) resonant plane waves (see eq. (HSJ)); this provides a set of amplitudes Aj, and finally 
one confirms that the superposition of the (Aj, ujj) plane waves fits the numerical solution 
*„(*)■ 

While immobile discrete breathers can be described as a sort of homoclinic (and time 
periodic) connection on the rest state, the mobile localized core instead connects a specific 
linear superposition of low amplitude nonlinear plane waves. One could say that the localized 
core needs for its motion to "surf over" a specific extended state of radiation: 

s-1 

(*backg)"(*) = A i exp ^ kn ~ ■ ( 35 ) 

3=0 

We note that among the members of the (s-parameter) continuous family of (p, q) reso- 
nant plane waves (see Section I), the fixed point solution contains only a particular member 
(Aj, ujj) from each branch (see Fig. ??.a). This selection varies smoothly with the (adia- 
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batic) continuation parameter v. In particular, the amplitude modulus \Aj\ selected increases 
smoothly from its zero value at the integrable limit (y — 0), for both signs of v. 
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FIG. 6: (a) Plot of the graphical solving of the resonant condition (in the Aj — > limit) for a (1, 2) 
resonant breather with Ub = 2.384 and Vf, = 0.189. (b) Power Spectrum S(co) of the background of 
this solution at v = 1.0. From (a) eq. (j!7j) gives the contribution of seven plane waves (j = 0, 6) 
but only five (j = 0, ...,4) of them are visible due to the difference of orders of magitude between 
the amplitudes \Aj\. The agreement between the resonant condition equation (for the fitted value 
of Aj) and the frequencies observed in S(ui) is up to machine accuracy. 

If the bare core of a fixed point solution (i.e. after substraction of the background) is 
taken as initial condition for a direct numerical integration of the equations of motion, one 
observes radiative losses, along with the corresponding changes in shape, velocity, etc. of the 
localized moving core. The motion of the bare localized core (not anymore a solution) excites 
extended states of the lattice. Thus, regarding the exact fixed point solution, one could say 
that radiative losses of the running core are exactly cancelled out when the localized core 
runs, with specific velocity, on top of the specific linear combination of (Aj, ujj) resonant 
plane waves (J33J). 

A complementary numerical observation is the following: Taking as initial condition for 
a direct integration of the equations of motion (J7|), a superposition of an immobile discrete 
breather and the background of a (p, q) resonant mobile breather, it evolves into a moving 
discrete breather, with approximate velocity Vb = (pub) / (27rq) . One thus would say that 
the background promotes breather translational motion with adequate velocity. In the 
next section, a connection between background characteristics and the particle perspective 
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(i.e. the Peierls-Nabarro barrier of collective variable theories), will be established, further 
illuminating the physical description of discrete breather mobility. 

Whatever physical perspective one may prefer, the numerical fact is that the generic 
structure of the fixed point solution is given by the superposition (J33j) . Not too far from 
v ~ 0, where the amplitudes Aj of the fixed point background have small values, one can 
carefully check that if the bare core is given as a starting guess for Newton iteration, this 
converges well to the exact complete solution (core + background), by developing the specific 
selection of Aj amplitudes. This confirms the robustness of the numerics. 

Though previous observations of nondecaying tails of numerically accurate mobile discrete 
breathers in Klein-Gordon lattices 20] and/or (solitary) traveling waves [31( in self-focusing 
equations had been reported (see also the interesting discussions on this issue in j2lj and 
|32j), no systematic study of those tails is known to us. However we clearly see that they 
are an essential part of the exact solution. As argued in the introductory section, the 
translational motion of a discrete breather introduces a new time scale. In a nonintegrable 
context, this fact unavoidably implies resonances with plane wave band spectra, and an 
exact self-sustained moving DB solution could only exist on top of a developed resonant 
background. This seems to have been (with a few exemptions) not fully appreciated in most 
of current litterature on mobile breathers, where the backgound is most often either ignored 
or deliberately suppressed. 

A notable feature of the plane wave content of the background $ Dac kg * s ^ na ^ ^ ne am ~ 
plitude modulus \Aj\ in differ by orders of magnitude, i.e. \Ai\ ^> \A 2 \ 3> | ^4_3 1 - - - , so 
that only a few frequencies are dominant for most practical purposes (see Fig. ??.b). In 
other words, the extended background associated to a spatially localized moving core is, in 
turn, strongly localized in the reciprocal (fc-space) lattice. The possible relevance of this 
observation is further discussed below in the concluding section. 

B. The background amplitude 

In order to characterize the specific features of the nonintegrable motion of discrete 
Schrodinger breathers, we focus here on the (perhaps) most remarkable among those fea- 
tures: the background amplitude of the uniquely continued fixed point. How does it evolve 
along the continuation path in parameter space? 
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For positive values of v we have followed the line in parameter space // + v = 1 (see 
equation (|7J)), while for negative values, we took the path ji — v = 1. We do not expect 
other paths to make important differences. As stated earlier, near v ~ 0, the amplitude 
grows from its zero value (at the integrable limit) for both signs of this parameter, for it 
is a nonintegrable effect. However, for larger values of nonintegrability \v\ the background 
amplitude evolution shows some important differences for the two signs of v. 

In Fig. [7|we plot the background amplitude (modulus) of the (1,1) resonant fixed point, 
versus the continuation parameter z/, for three different values of the breather frequency oui,. 
For v > 0, one observes that the amplitude steadily increases with v before continuation 
stops (i.e. Newton iteration ceases to converge beyond a certain maximum v value). Note 
that the amplitude grows faster for higher values of the frequency, and that the continuation 
stops (correspondingly) at a smaller value of v. This may suggest that the failure of fixed 
point continuation is related to a somewhat excessive growth of the background amplitude, 
an issue that will be discussed later. 

For v < 0, after an initial growth the background amplitude decreases down to almost 
negligible values around v ~ —0.3, then grows and again decreases close to zero at v ~ —0.39, 
and so on, in progressively narrower intervals with larger peak amplitude, until continuation 
stops. Most noticeable is the fact that the intervals neither depend on the breather frequency 
u b nor on the breather velocity v^. Why do background amplitudes decay so dramatically 
at those regions in parameter space? An important hint is presented in the next section, 
where the Floquet stability analysis of immobile discrete breathers will show a coincident 
situation of mirror-symmetry breaking (and its absence for positive v values). 

For other values of p and q that we have numerically investigated, the same features of 
the background amplitude variation as shown in Fig. [7| are qualitatively reproduced. 



C. Floquet analysis 



On the basis of the general arguments given in [27, 33], the Floquet spectra of immobile 
DB in the thermodynamic limit, N —>■ oo, consists of two components: the (continuous) 
Floquet spectrum of the asymptotic state of the solution (rest state), and a discrete part 
associated with spatially localized eigenvectors. The continuous part is composed by small 
amplitude (linear) plane waves, the so-called phonons. However, for mobile DB the asymp- 
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FIG. 7: Background amplitude versus v for three different (1/1) resonant breathers with frequen- 
cies: (a) ujb = 5.65, (b) u>b = 4.91, (c) u>i, = 4.34. Note the two different behaviours: for positive 
values of v \<&backg\ 2 is a monotonous increasing function of v while for the negative part it shows 
smooth rises and falls. 

totic state of a (p, q) resonant fixed point solution is a superposition of plane waves, the 
background $backg- From this, one should expect the Floquet spectrum of a (p, q) resonant 
DB being composed of two components: the discrete (spatially localized eigenvectors) and 
a continuous part associated with the linear stability of the background plane waves. The 
continuous part of the Floquet spectrum should reflect the same results of the modulational 
instabilty analysis of section Hi CI In particular, this means that any modulational instability 
a plane wave may suffer will be also an instability of a fixed point solution whose background 
contains this plane wave. In the future we will refer to to any instability of the continuous 
part of the Floquet spectrum as background instability. Any instability from the discrete 
part is a core instability. 

First we focus on core instabilities. For this we turn attention to the continuation of 
mobile (p, q) resonant breathers. Fig. |H] shows in the v — u b plane (dotted line), the values 
Vmaxi^b) where the numerical continuations stop due to nonconvergence of Newton iteration 
for p — 1, q — 1 and v > 0. As was remarked above, the continuation stop is associated 
with the rapid increase of the background amplitude shown in Fig. El Only low frequency 
breathers, for which the background amplitude increases more slowly, can be numerically 
continued all the way to the standard DNLS equation. The linear stability analysis of (p, q) 
resonant breathers yields a well defined region in the v — u>b diagram where core instabilities 
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appear. There is an island inside the continuation region of Fig. |H1 where the Floquet 
spectra contain a real eigenvalue A > 1. We observe the evolution of this Floquet eigenvalue 
(and its complex conjugate) as the parameter v is increased in Fig. 0a, for a (1, 1) breather 
of frequency u>b = 2.678. Here the angle {Opioq) in the complex plane is plotted versus v. 
The interval of constant zero angle corresponds to the section (constant u>b) of the instability 
island in Fig. |H| 




0.2 0.4 0.6 0.8 1 



v 

FIG. 8: Continuation diagram of (1, 1) resonant breathers as a functon of the frequency uj^- The 
end of the numerical continuation, fmaxi^b) , is represented by the line with dots. The region where 
mobile breathers suffer from core instabilities is limited by the shaded area. 

Along the whole continuation path, the profile of the corresponding unstable eigenvector 
is localized. An example of this profile inside the instability island is shown in Figs. Ob and 
Elc, where one observes that the localized instability shows a decaying background along the 
direction opposite to the motion. The decay rate increases as the modulus of the eigenvalue 
grows and decreases again when A returns to the unit circle. On the other hand, the stable 
Floquet eigenvector associated with 1/A shows a wing decaying along the mirror symmetric 
direction. The direct integration of the equation of motion reveals that the unstable solution 
experiences a pinning after a transient of regular motion with velocity Vf, = p/{qTb). After 
the solution pins at site n, its core center oscillates around this site. The trapping of the 
unstable MB could be interpreted as a result of the energy losses that the growth of the 
linearly unstable perturbation induces on the solution. 

Returning to the instability island shown in the diagram of Fig. |Hl some final observa- 
tions are worth summarizing: (i) there is a range of frequencies where mobile breathers of 
the standard DNLS equation [y = 1) suffer from this instability; (m) very high frequency 
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breathers do not experience this instability (in the short range where they can be continued); 
(Hi) very low frequency breathers are stable all the way up to v — 1. 
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FIG. 9: (a) Floquet angle evolution of the spectra of a (1, 1) resonant breather with lu^ = 2.678. 
The thick trajectory corresponds to the localized eigenvector that becomes unstable (0Fioq = 
interval). Instantaneous profile of the real (b) and imaginary (c) part of the Floquet unstable 
eigenvector of a (1, 1) breather with ujb = 3.207 and v = 0.26. The decaying tails along the 
direction opposite to the motion reveals the energy loss that the unstable eigenvector causes to the 
solution. 

We turn now to background instabilities. Once we know the plane wave content (A; , 
ki,..) of a (p/g)-resonant fixed point, we can know whether the solution is subject to MI 
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or not and, if it is unstable, what are the harmful perturbations (Q). This problem is not 
so simple because we cannot know a priori the plane wave content if we do not have the 
amplitudes of each one (|17|) . However, we can derive a necessary condition for not having 
MI if we consider that, from (J17|) . the background is always composed of at least one plane 
wave (m = 0) with k between [— n/2,0]. From this we can simplify the analysis of the 
background stability to the ko plane wave stability as a necessary condition for the MB 
stability For this we calculate, for each v and k, the value of the right-hand side of (f2*rjj) for 
all the range of Q ( [— tt, 7r] ) and A. If this value is always positive the plane wave with this 
ko is free from modulational instabilities at this point of the model (JJJ) with parameter v. 
From this extensive exploration we obtain, see Fig. EH the region in the k — v plane where 
MI is present. 
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FIG. 10: Modulational instability existence diagram for a plane wave with wave number ko £ 
[— 7r/2,0]. This diagram fixes the region where mobile discrete breathers with a background com- 
posed of only one plane wave do not suffer from background instability. 

In the range of v between [—1, —0.5] there is no modulational instability for single plane 
waves of any value of k between [— 7r/2,0], and in particular for k . However, this does 
not guarantee that moving breathers are free from these instabilities in this region, unless 
the background has only one plane wave (as is sometimes the case). On the contrary, in 
the region v > any moving breather suffers such instabilities. The transition area in the 
region v G [—0.5, 0] presents MI depending on which ko we have. For the range where no 
plane- wave with k between [— 7r/2, 0] is subject to MI we can assure that if there is only one 
contribution, k , to the background the corresponding MB solution is stable. For example, 
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Unstable 
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this is the case for (1/1) resonant breathers if ^ > 4 and for (1/2) resonant breathers if 
ujb > 8.46. The Floquet spectra of a moving breather satisfying these requirements is plotted 
in Fig. ITTlc. 



(a) (b) (c) 
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FIG. 11: Floquet spectra of (1, 1) resonant breathers: (a) for ujf, = 4.348, Vf, = 0.692 and v = 0.08 
the spectra shows the core (localized) instability; (b) for uj^ = 6.610, = 1.052 and v = 0.07 the 
spectra shows the background (modulational) instability (also present but not visible in (a)); (c) 
for u>b = 4.348, = 0.692 and v = —0.39 the solution is linearly stable. 

After the analysis of both types of instabilities eventually experienced by moving 
Schrodinger breathers, we finally report on a most relevant numerical fact revealed by the 
Floquet analysis of the family of immobile discrete breathers for v < 0: Near v ~ —0.3 
an immobile two-site DB experiences a mirror symmetry-breaking (pitchfork) bifurcation 
becoming linearly unstable. When approaching the bifurcation point, two conjugate Flo- 
quet eigenvalues quickly approach +1, where they meet, and then separate along the real 
axis. The eigenvector associated to the unstable A > 1 Floquet eigenvalue is localized and 
odd-symmetric, and is termed the symmetry-breaking or depinning mode <j) dep . We recall 
here that the background of an immobile breather is the rest state $ = 0, whose continuous 
spectrum consists of small amplitude (linear) plane waves. The depinning mode, on the 
other hand, is a localized core instability of the immobile breather, favoring a translation of 
the core center. For a smaller value of v ~ —0.39 there is another symmetry-breaking bifur- 
cation where the two-site breather becomes stable, again interchanging the stable character 
with the one-site. The corresponding bifurcation diagram for these two symmetry breaking 
transitions is plotted in Fig. [121 

In the first symmetry breaking bifurcation, two unstable mirror-asymmetric immobile 
breathers emerge from the bifurcation point, progressively evolve toward the (stable) two- 
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FIG. 12: Graphical representation of the two first symmetry breaking bifurcactions for v < 0. 
The quantity £ in the vertical axes of both figures is defined, referred to the one-site breather, as 
the difference between the modulus |$| of the two sites adjacent to the maximum (l^mazl), i-e. 
£ ~ l^maz-il — I'&maaj+il- For the one-site DB ^ = and for the two-site DB £ = 1, for this £ is 
conveniently normalized with the the difference between $ max and $ ma x±i- The continuous lines 
represent the regions where the static solutions are linearly stable while the discontinuous ones 
represent the unstable regions. The modulus profile of the three immobile coexisting solutions are 
plotted in the central insets for uib = 6.215 and v = —0.3012. 



site breather, and finally collide in a new pitchfork bifurcation from where a unstable two-site 
breather emerges. The net result is an inversion of stability between one-site and two-site 
immobile breathers. Around the narrow interval of v values where these two bifurcations 
occur, the energies of the three types of breathers involved (one-site, two-site, and asymmet- 
ric) have very small differences. From a particle perspective, this should make the breather 
motion easier. It is precisely in this same narrow interval where (see IIIIB|) we observe that 
the background amplitude of moving breathers becomes negligible. This is not a coincidence 
as we will argue in II VI 



IV. PARTICLE PERSPECTIVE ON DISCRETE BREATHERS 



The appealing framework and success of collective variable approaches (see e.g. 



y,y, 



HQ, 



to the problem of nonintegrable motion of discrete breathers relyes on the fidelity 
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of a particle-like description of these field excitations that they provide. In these approaches, 
the effective dynamics of only a few degrees of freedom (e.g. the localization center, and the 



12j ) replaces the whole description of the 



spatial width of the state, etc... in some instances 
moving localized state. 

Though unable to account for all the nonintegrable features, perturbative collective vari- 
able theories of NLS lattices provide a sensible physical characterization of important fea- 
tures of the nonintegrable mobility of localized solutions, like the emergence 3^ of a Peierls- 
Nabarro barrier to motion. Here we summarize the main results of this particle-like descrip- 
tion and compare them with the behaviour of numerically exact (p, q) resonant moving 
breathers. Our goal is twofold: to acquire a correct physical understanding of the numer- 
ical facts, and then to make an assesment of validity and intrinsic limitations of collective 
variable approaches. 



A. Collective variables theory. 



A presentation of the particle perspective on moving Schrodinger breathers near the A-L 
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where the interested reader 



integrable limit can be found in 36J (see also 
will find the relevant formal aspects of the theory. 

Using the integrable solitary wave (jlj) as an ansatz for the moving breather solution in 
the perturbed A-L lattice, v ^ and small in (|ZJ), one considers the parameters a, /3, xq and 
Q as dynamical variables (variation of constants). The time evolution of these parameters 
in the perturbed lattice is governed by: 

sinh f3 



where 



x 
a 



2 sin a- 



P 



2 cos a cosh (3 + aio + g{(3) 


8tt 3 sinh 2 (3 



—v 



(3 3 sinh(7r 2 s//?) 



sin(27rsa;o) 



(36) 
(37) 
(38) 

(39) 



r 2 sinh (3 cosh (3 sinh 2 (3 
2u\ ^ - 1] 



+ V 4-7T 2 COs(27TgXo) [ 



s=l 



sinh 2 (3 cosh(7r 2 s / (3)tt 2 s 
/3 4 sinh 2 (7r 2 s//5) 



27 



2 sinh 2 /3 2 sinh (3 cosh (3 



p 3 smh(7r 2 s/(3) (3 2 smh(ir 2 s/f3)" 

These relations can be viewed as the Euler-Lagrange equations of the collective variable 

I I 

Lagrangian obtained in |36[. The variation of the breather parameters give the evolution of 
solution (|3]) for the perturbed A-L equation. Furthermore, one can regard eqs. (fHHj) and (jSHJ), 
as the Hamilton equations for the canonical conjugate variables xq and a of the following 
effective hamiltonian: 

_ sinh (3 

rtcS = Teff + Veff = 2 COS OL — — 

^ 47r 2 sinh 2 /? 

2^ A3 • u? 2 7m cos 2ttsx , 41 
^— j' p d smh^s/p) 
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This effective hamiltonian dictates the dynamics of the position of the solitary wave. 
Note that the (collective) variable f3 is an invariant of motion, so it enters as a parameter 
into the effective hamiltonian, and that the time-average value of Cl (the parameter uj^ of the 
integrable solitary wave, now a function of time) is an increasing function of this parameter 
(3. The effective potential V e fr acts as a barrier to the displacement motion (x variations) 
and is naturally related to the Peierls-Nabarro potential. The amplitude of this barrier is 
an increasing function of both the nonintegrability parameter \v\ and (3. The equilibrium 
points (representing immobile breathers) of this potential are Xq = n and n ± 1/2 with n an 
integer. For a = 0, the former are stable (centers) one-site breathers, while the latter are 
unstable (saddle) two-site breathers; for the case a = n (staggered breathers) the stability 
is reversed. 

A remarkable further consequence is the following 35]: there are no perturbative traveling 
wave solutions, for values of v larger than certain critical value v C r{!3)- In particular, for 
/3 > /3 C ~ 3.6862 , one cannot continue A-L mobile breathers (i.e. v CT = (see also important 
remarks in 



[38]). This consequence could be also (qualitatively) expected for a class of 
nonintegrable Schrodinger lattices (for some qualified perturbations of the integrable limit) 
with on-site nonlinearity. One expects also that lattices with purely inter-site (FPU-like) 
nonlinearity do not show this kind of transition. 

In Fig. El we plot typical phase portraits at both sides of v cr . Fig. El a shows the 
dynamics for v smaller than the treshold value (given by (3): there are open trajectories 
in Xq corresponding to mobile breathers and closed orbits between the separatrix manifolds 
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corresponding to breathers which oscillate around the equilibrium position of V e s- Fig. [T31 b 
is the phase portrait after the transition: there are no longer mobile solutions and (besides 
the oscillating breathers) there are instead open trajectories in a. The transition point, for a 
given (3, occurs when trajectories with rotating a appear, and moving breathers disappear as 
the effect of separatrix line rearrangement on the cylinder (xq, a(modulo27r)) phase portrait. 

(a) (b) 
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FIG. 13: Collective variable (a, xq) phase portrait transition for a value of /3 = 3.0. (a) Shows 
the phase portrait for v = 0.2 (< v cr \ there are xo-unbounded trajectories (mobile breathers) 
coexisting with bounded ones (oscillating breathers). When v = 1.0 (> v cr ) (b) we only have 
xo-bounded trajectories: there are no mobile solutions. 

Note that the existence of oscillating breathers is a consequence of the existence of a 
Peierls-Nabarro potential. These breather solutions do not perturbatively continue from the 
integrable limit. In II V CI we will investigate them and provide further numerical confirmation 
of the existence of these genuinely nonperturbative solutions, predicted by the collective 
variables theory. 

B. Energy balance governs mobility. 

In order to correlate collective variable predictions with the numerical results presented 
in section HID one should first realize that our direct numerical approach computes breathers 
with fixed values of Ub and f& and that these parameters are not tied to any specific ansatz. 
In particular, the connection of these two parameters with the collective variables is given 
by eq. (jHJ) in the integrable limit. For the perturbed (near-integrable) lattice, ujf, and Vb are 
identified as the time averages of fl and xo, respectively. 
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The Peierls-Nabarro (PN for short) barrier is naturally identified as the energy difference 
(given by the Hamiltonian (jHJ)) between the two immobile breathers of the same frequency 
u b , one centered at a site n and the other (two-site) at a bond n ± 1/2 : 

E PN (u, u b ) = H(u, u; b , n) - H(u, u b , n±^) (42) 

In the integrable A-L limit this barrier is zero due to the degeneracy (continuous translation 
invariance) of the breather family solution, but for v ^ this invariance is broken and only 
these two isolated solutions persist. The energy difference of the two pinned solutions is thus 
viewed as the minimal extra "kinetic energy " of center of mass translation that a mobile 
breather must incorporate for overcoming the barriers to its motion. 
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FIG. 14: Density Plot of the absolute value of the Peierls Nabarro barrier, \Epn\, as a function 
of ui b and v. For positive values of v, \Epn\ is a monotonous increasing function of v and u. For 
negative values the plot reveals the oscillating behaviour of \Ep^\ as a function of v (for a given 
value of uj b ). 

We have studied the behavior of the PN barrier in the Salerno model by continuing 
immobile breathers, both centered at a site and at a bond, while computing their energy 
difference. The computations of the barrier are made for a grid of values of u b . Fig. El 
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shows the " equipotentiaT lines of the PN barrier in the (is, Ub) plane. The results show 
different behaviors depending on the sign of is: 

is < 0. Here one observes the effects of the symmetry-breaking bifurcations cascade de- 
scribed in 1111 CI The successive stability inversions between site and bond centered breathers 
involve a substantial decrease of the Peierls barrier. The appearance of asymmetric solu- 
tions in these bifurcations introduces a new energy and, correspondingly, the Peierls barrier 
is computed as the maximum energy difference between the three pinned solutions: the two 
symmetric (site and bond centered) and the asymmetric breather. 

is > 0. Here the behavior of the Peierls barrier follows qualitatively the collective variable 
predictions on the effective potential experienced by the particle. The increasing character, 
with is and W&, of the numerical barrier is qualitatively the same as that predicted from V e ff 
(as a function of is and j3) by the theory. 

The PN barrier of (u;&) immobile breathers and the background amplitude of (oj^, Vb = 
mobile breathers are in fact strongly correlated. This correlation is obtained considering the 
functions \Epn\(is) and \&backg\ 2 (v) ■ Both functions are plotted for a fixed value of uJb = 4.34 
in Fig. E3a. The behavior of \Epn\(is) for negative is (revealing the cascade of bifurcations 
explained before in IIII C|) is closely followed by \&backg\ 2 (v) with the corresponding sequence 
of growths and decays. The strong correlation holds also for positive values of is, where 
numerical PN barrier data are available for a larger interval of v values (due to the absence 
of the symmetry-breaking cascade of bifurcations). Indeed, the correlation is so strong that 
one is tempted to view the PN barrier and the background amplitude as complementary 
aspects of a single phenomenon: the breaking of the continuous translational invariance, and 
the associated lack of core momentum conservation (^(J. Indeed, the background amplitude 
of moving breathers is a monotone increasing function of the PN barrier of pinned breathers 
of the same frequency, as shown in Fig. H51b. where \&back g \ 2 (\EpN I) is plotted. 

However, we also observe clearly in Fig. ITKI a that, when the continuation end is ap- 
proached, the rate of growth of \&back g \ 2 (v) increases dramatically (the concavity of the 
curve in log scale turns upwards), while the PN barrier does not increase much faster than 
before. This is reflected in figure fI31b. where the slope approaches verticality, indicating 
that, in this range of Epjy values, the background grows rapidly. 

This numerical observation suggests taking a closer look at the precise influence of the 
background amplitude on the core energy variations associated with the existence of PN 
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FIG. 15: Peierls Nabarro barier |-Epjv| from inmobile breather with cj& = 4.34 and background 
square amplitude \&backg\ 2 for a (1,1) resonant breather of the same frequency (vb = 0.691). In 
(a) we show both quantities in semi-log scale as functions of v, illustrating the strong correlation 
between them for both signs of v. Figure (b) shows for positive values of u, the nonlinear relation 
between |$6 ac A; g | 2 and \Epjy\. Note the sudden increase of the slope close to the end of numerical 
continuation. 

barriers. To this end, we use the conservation of the Hamiltonian (JHJ) and insert this equation 
into the form (jHHj) of the (p, q) resonant fixed point. The energy of the solution can be 
decomposed in the following terms: 

n[^] = Hficare] + «[l^ 9 ] + H mt , (43) 

where TC mt is the interaction energy, i.e. the crossed terms of & COTe and &back g - Let us now 
consider the simplest case in which the background has a single resonant plane wave. Along 
with the total energy, also the energy of the plane wave is a constant in time so that 

-^r~ = w ■ (44) 

In other words, the variations of the core energy during the motion are exactly compensated 
by those of the interaction term. 

If one takes, as an ansatz for $ core , the A-L solution, one formally obtains for Ji COTe = 
7~i\§core\ the collective variables Hamiltonian (|41jl . But note that here it would not any- 
more be a constant of motion, due to the interaction with the background. Instead (as 
P. Kevrekidis suggested to us), we directly compute numerically the evolution of the core 
energy 7i core , which in turn determines 7i mt up to an additive constant. 
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For this we take a fixed point solution with a single plane wave in its background, and 
then substract off the plane wave to obtain $ cor . e , from where TC core (t) is computed. In Fig. 
[Tolwe have plotted the evolution of the core energy as a function of the localization position 
(center) of the breather core. The localization center of a lattice function $ n is defined using 
the conserved norm (JDJ): 

_ E n nln(l + /i|$ n | 2 ) 

JJ7 • (45) 

As expected, the core has extracted the maximum available from the interaction energy 
(with the background) when the core passes at n ± 1/2 (maxima of the PN barrier) and has 
returned it to the interaction term when at n (minima of the PN barrier). 
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FIG. 16: Plot of H core of a (1, 1) resonant breather as a function of the localization center xq for 
different values of v. The parameter of the solution are uji, = 5.056 and Vb = 0.805. The values of 
v are 0.04, 0.08, 0.12, 0.16, 0.20, 0.24, 0.25 and 0.2512 (end of the continuation), the amplitude of 
the oscillation of 7i mt grows with v. (The minimum value of 7i tnt has been set to zero in order to 
compare the differents functions.) 

Another interesting feature of these numerically obtained functions is seen from the vari- 
ations in the form of the oscillation of 'H core as the nonintegrable parameter v is increased. 
At the same time, as the energy difference between n and n ± 1/2 increases the modulus of 
the derivative d7i core / dx$ in the neighborhood of Xq = n also increases. These variations 
become faster when the end of the continuation is approached, reaching a cuspidal point for 
the last v reached. The background amplitude is included in H mt , and of course in &H mt /dt; 
the dramatic variation of it at the end of the continuation could be interpreted in terms of 
this derivative variation in xq = n. 
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C. Oscillating breathers 



The emergence of the Peierls barrier and the behavior of the background amplitude 
illustrate the physical interpretation of this background as a (p/q)-resonant energy support 
to overcome the barrier to motion. We now confirm this statement searching for another 
kind of solution: oscillating breathers. These solutions are predicted by collective coordinates 
approaches and are a consequence of the loss of translational invariance out of the integrable 
limit. Following the above interpretation of the background role one can imagine certain 
solutions with a background amplitude not high enough for surpassing the Peierls barrier 
and allowing travel along the lattice. In terms of a well defined potential, considering the 
particle perspective, the center of these localised solutions would be oscillating between 
(n — 1/2) and (n + 1/2) for the unstaggered ones or between n and (n ± 1) for the staggered 
ones. 

From our perspective, the oscillating breathers are solutions with two frequencies: the 
internal one of the breather (u;&) and the one corresponding to the oscillatory motion (u> sc) ■ 
Once again, we have a problem dealing with two time scales and consequently we have to 
impose that the two frequencies are commensurate pu b = qu osc . The fixed point problem is 
now associated with the map: 

T qTb $ n (t) = $ n (t) (46) 

We cannot, however, develop the Newton iteration scheme in a similar way as for mobile 
breathers. There is no longer any family of oscillating breathers providing a good start point 
for the continuation (they are intrinsic solutions of the nonintegrable regime because they 
appear as the Peierls barrier emerges). The way to obtain a good ansantz (as Cretegny and 
Aubry already used to find mobile breathers in Klein-Gordon lattices j^J) is to perform a 
small perturbation of the static solution (pinned at a site n) with the depinning internal 
mode: 

^ansantz = ^static ^ + ^dep ^ 

The dynamics of the perturbed solution for small enough values of e shows the oscillating 
behavior expected and for large enough values of e the breather starts to move. Obviously 
in both cases the motion finishes after a transient due to radiation (they are not exact 
solutions). Tuning the parameter e we search for those oscillatory transients whose u osc 
is resonant with the breather frequency u^. The transient is much more stable when the 
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nonintegrable parameter v is very small, close to the A-L limit. We first search here for a 
good initial guess for the method and then obtain the exact solution of the map ([46)1 . Once 
the exact solution is obtained for a small u, we can perform the continuation to higher values 
in the same way as we did for mobile solutions. In Fig. \n\a, we show the evolution of the 
amplitude of oscillation as v is increased from 0.05 to 0.18. The amplitude of the oscillation 
is represented by the phase portrait of the localization center of the breather defined as in 
f!45|) . The continuation reflects that the amplitude of the oscillation, for a fixed value of u osc , 
grows with v. In Fig. Elb the density plot of |$ n | 2 is shown as a function of time, revealing 
the oscillating pattern of the solution. 

The existence of exact oscillating breathers is a consequence of the existence of a Peierls 
barrier. The structure of these solutions reveals the existence of a background (resonant 
with the map) whose amplitude grows as v (and consequently the amplitude of oscilla- 
tion) is increased. This is the picture we expected from the role played by the interaction 
background-core in the energy balance during motion. The monotonous growing behavior of 
the background versus the oscillation amplitude, strongly suggests that if the amplitude of 
the former is increased the solution will be able to translate steadily. This has been checked 
by direct numerical integration, because no exact solutions connecting the oscillating with 
the mobile ones can be obtained due to the different maps employed to obtain both types of 
solutions. However, the existence of a background in the exact oscillating breather solutions 
and its behavior with the amplitude of the breather oscillations are fully consistent with the 
interpretation of the results obtained for the mobile solutions. 



D. Validity and limitations of particle perspective 

The most basic result of the perturbative collective variable theories away from the inte- 
grable regime is the existence of a Peierls-Nabarro potential function of the core (collective 
variable) center. It expresses (in particle-like terms) that the breather position is no longer 
indifferent because the continuous translational invariance has been broken. From this also 
naturally comes the existence of oscillating breathers. We have seen how our numerics fully 
confirm the qualitative validity of these predictions. 

A further prediction concerns the phase portrait's transition studied in j^. Despite 
the fact that our end of continuation is correlated with the equipotential lines profile of 
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FIG. 17: (a) Evolution of the localization center xq of an exact (l/18)-oscillating breather for 
different values of v: 0.05, 0.06,.., 0.18. The internal frequency is u>6 = 3.086. The amplitude of the 
oscillation of xq increases with v revealing the nonlinear character of the motion for the highest 
values of v. (b) Density plot of the time evolution of |$ n | 2 for the above oscillating breather. 

the numerical PN barriers, and the phase portrait transition is also related to their sudden 
growth, no clear connection (between transition and end of continuation) can be established. 
The end of continuation is itself sensibly interpreted as a numerical consequence of the 
sudden increases of the amplitude background, and does not imply neccesarily the existence 
of any global phase portrait transition. 

However, in some respects the perturbative collective variable theory is clearly incomplete: 
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For example, it is unable to predict the observed localized (core) instability bifurcation and 
the observed symmetry breaking transitions for v < 0. These bifurcations could easily 
appear in a theory with (at least) two variables (a dimer) experiencing the Peierls-Nabarro 
potential, which would demand an improved perturbative ansatz. This improved ansatz 
must coincide in the integrable limit with the A-L solution. One can use the numerical 
results to guide the construction of such an improved ansatz. In this respect, the following 
observation may be relevant. The parameter f3 of the A-L solution determines both the 
amplitude and the width of the localized pulse. However, our numerical estimates of these 
breather characteristics for immobile breathers show clearly that, for fixed value of uib, the 
breather width is independent of v, while the amplitude varies with it. In other words, away 
from integrability, width and amplitude of the (immobile) breather are no longer a single 
collective variable. 

Beyond any other limitation of the perturbative collective variable theory, the back- 
ground (an indispensable part of the exact solution) is absent in the perturbative ansatz, 
and it cannot appear later in that context. A complete theory of (nonlinear Schrodinger) 
breather motion should somehow incorporate the background in the ansatz itself. If correct, 
it should then predict that the background amplitude grows from zero with the nonintegra- 
bility parameter z/, and (ideally) so on with all the numerically observed behaviors. One 
possible way to develop the analytical approach could be to use the method presented in 
4l| . In this scheme, eq. (jB}) may play an important role, for it provides the energy balance 
governing the translational motion of the breather core. In other words, our results show 
that the core energy is not an invariant of motion and this requires the existence of a finely 
tuned background whose nonlinear interaction with the core compensate the core energy 
variations. We hope that the numerical work presented here will provide motivation for 
further analytical research. 

V. CONCLUSIONS AND PROSPECTIVE REMARKS 

We have studied numerically the problem of mobility of nonlinear localized solutions of 
the one-dimensional Nonlinear Schrodinger lattice ((Zj) using a regularized Newton method. 
This method allows us to continue the family of mobile Ablowitz-Ladik discrete breathers 
into the nonintegrable domain of model parameters. We find that these solutions decay 
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asymptotically, in space, to an excited lattice extended state (the background), whose am- 
plitude vanishes at the integrable Ablowitz-Ladik limit. This component of the solution is 
unambiguously found to be a linear combination of nonlinear resonant plane waves whose 
amplitudes decay typically, in /c-space, exponentially. The exponentially localized oscillation 
(the core) of the amplitude probability rides over this extended radiation state. 

$ = $ core + $ backg (48) 

Numerically exact moving discrete breathers with an infinitely extended tail of small am- 
plitude were already observed in some cases for Klein-Gordon lattices with Morse potential 
by Cretegny and Aubry j^J, however no investigation of the background of these exact 
solutions is reported, so they were able only to "..suggest that generally a strictly localized 
breather cannot propagate without radiating energy". Our systematic study of the NLS 
lattice allows us to go further by showing that the extended background (here fully charac- 
terized) plays an important and subtle role in the translational motion of the localized core. 
Indeed, it is an indispensable part of the exact solution in the nonintegrable regime. Ex- 
act mobile localization only exists over finely tuned extended states of the nonlinear lattice. 
Mobile "pure" (i.e. rest state background) localization must be regarded as very exceptional 

The high numerical accuracy of the method allows a precise Floquet stability analysis of 
the moving discrete breather solutions. The Floquet analysis reveals the type of instabilities 
(both localized and extended) eventually experienced by the two-parameter family of moving 
breathers. Some generic bifurcations in the space of model parameters are thus identified. 
For negative values of the nonintegrability parameter u, we find narrow regions where the 
immobile breathers experience mirror symmetry-breaking bifurcations and, simultaneously, 
the amplitude of the background component of the mobile breather solutions decreases down 
to almost negligible values. 

The most relevant predictions of perturvative collective variable theory are confirmed 
by our numerical results, which show the existence of Peierls-Nabarro barriers to breather 
translational motion. Furthermore, the existence of exact oscillating breather solutions is 
numerically confirmed. They are found to contain an extended background whose amplitude 
is typically much smaller than for mobile breathers. 

The Peierls-Nabarro barrier Ep^, computed from immobile breathers, and the amplitude 

38 



background of moving breathers are found to be strongly correlated, which correctly sug- 
gests that the background has a role in the energy balance required to overcome the barriers 
to translational motion. This is also fully consistent with the observations on the back- 
ground amplitude behaviour of spatially oscillating anchored breathers. Indeed, assuming 
the simplest case of a monochromatic background, the variation of the core energy during 
its translational motion is exactly compensated by the variation of the energy of interaction 
between the localized core and the background. Currently used effective particle (collective 
variable) theories are thus seen as intrinsically incomplete, because core energy is not an 
invariant of motion. Any sensible improved approach must adopt equation (|48jl as starting 
point for improved perturbative ansatze, and we hope that our work will estimulate further 
studies along these lines. 

There are, at very different levels, several open questions to further research. From a 
technical point of view, it is important to analyze carefully the irrational limit p/q — > a, 
of the solutions. In particular, in this limit the number of resonant plane wave branches 
tends to a continuum and one could (or not) expect that exponential localization in the 
reciprocal lattice persists in that limit. This can be addressed numerically, though systematic 
investigations may require some efforts in optimizing the time efficiency of current numerical 
schemes. 

An important issue regarding applications is the phenomenology of multibreather states. 
In particular, studies on collisions of a pair of breathers may find in this study of exact mo- 
bility a useful reference in order to deal with the complexities that emerge from the many 
time-length scales involved in these physically relevant phenomena. Much simpler multi- 
breather states, e.g. train-like chains of (moderately) separated moving breathers could also 
be investigated. Not least, the perspective and results presented here may be of some interest 
to studies of the effects of coupling to (nonthermal and/or thermal) radiation baths in the 
breather and multibreather states of nonlinear lattices |42[ and the practical manipulation 
and patterning of localized "hot spots" by external fields [43j . 
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